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Antiferromagnetic condensates are generally believed not to display modulational instability and 
subsequent spin-domain formation. Here we demonstrate that in the presence of a homogeneous 
magnetic field antiferromagnetic spin-1 Bose-Einstein condensates can undergo spatial modulational 
instability followed by the subsequent generation of spin domains. Employing numerical simulations 
for realistic conditions, we show how this novel effect can be observed in sodium condensates confined 
in an optical trap. Finally, we link this instability and spin-domain formation with stationary modes 
of the condensate. 
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I. INTRODUCTION 

The appearance of spin degrees of freedom in atomic 
matter waves opens up possibilities for new phenomena 
such as spin waves [l[ , spontaneous magnetization and 
spin mixing Q. However, perhaps the most intriguing ef- 
fect is associated with_complex patterns, such as spin tex- 
tures 0] or domains ||, which may appear either as sta- 
tionary low-energy states or emerge spontaneously due to 
condensate instabilities. Pattern formation is a common 
feature in the dynamics of extended nonlinear systems 
ranging from optics p to fluids 7]. Such patterns of- 
ten develop through the exponential growth of unstable 
spatial modulations, known as modulational instability. 
In the spinor condensates we have the opportunity to 
examine such effects in an environment which is remark- 
ably easy to control and manipulate, simply through the 
addition of an external magnetic field. 

The origin of the intriguing physics of spinor conden- 
sates lies in the spin interaction between atoms, which 
allows for an exchange of atoms between different spin 
components. The parametric nature of this interaction 
mirrors similar effects observed in nonlinear optics, where 
the interaction of several optical modes may lead to the 
development of new frequencies jfj}]. Of particular inter- 
est to our case is the possibility that instabilities of an 
intense light beam may occur even when the wave is cou- 
pled to a spatially stable eigenmode and propagates in 
the normal-dispersion regime 8]; in this case the inter- 
play of natural and self-induced birefringence leads to 
nonlinear polarization symmetry breaking and polariza- 
tion modulational instability. By analogy we thus might 
expect similar instabilities in an initially stable polar con- 
densate subjected to additional spin component coupling 
through an external magnetic field. 

In the absence of an external magnetic field, the de- 
velopment of spatial modulational instability in three- 
component (or spin-1) ferromagnetic condensates and the 
subsequent formation of spin domains has been well es- 
tablished both theoretically [§, [l(J [H| and experimen- 
tally [Hj]. However, early work on the zero field case @ 
determined that antiferromagnetic (or polar) conden- 



sates are modulationally stable. Experimental observa- 
tions suggested this is also true for a weak magnetic field, 
however these experiments were carried out with a con- 
densate smaller than a spin domain [l3j . 

In this paper we reveal that in fact the presence of 
a weak magnetic field (~ 175 mG) leads to spin domain 
formation in antiferromagnetic condensates, provided the 
condensate is larger than the spin healing length. Fur- 
thermore we show that this spin domain formation is 
initiated by a new type of modulational instability, rem- 
iniscent of instabilities observed in nonlinear optics Q 
and not seen before in Bose-Einstein condensates. While 
spin-domain formation in antiferromagnetic condensates 
has been observed before in the presence of a magnetic 
field gradient Q, we show here that it occurs equally 
well in the presence of a homogeneous magnetic field. 
Furthermore we reveal that this modulational instabil- 
ity and spontaneous spin-domain formation is associated 
with stationary states which exist in the presence of the 
weak magnetic field, and which intrinsically break the va- 
lidity of the single- mode approximation (as seen earlier 
in [1J|). We discuss realistic experimental conditions for 
the observation of these novel effects. 

The paper is organized as follows. SectionHTlintroduces 
theoretical model of spin-1 condensate in a homogenous 
magnetic field. In Sec. lIIII we investigate homogenous sta- 
tionary states in magnetic field and analyze their stability 
with respect to plane wave perturbations (modulational 
stability). Section [TVl presents results of numerical simu- 
lations corresponding to experimentally relevant conden- 
sate evolution, demonstrating the possibility of observa- 
tion of new instability in antiferromagnetic condesate. In 
Sec. [V] we fink this instability and spin-domain formation 
with stationary modes of the condensate, and Sec. IVII 
conludes the paper. 



II. MODEL 

The evolution of a dilute spin-1 (F = 1) Bose-Einstein 
condensate (BEC) in a homogeneous magnetic field is 
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given by the coupled Gross-Pitaevskii equations 
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n = n+ + no + n._, and V(r) is an external potential. 
The nonlinear coefficients are: cq = 4tt?i (2a 2 + o )/3 
and c 2 — 47r?i 2 (a 2 — aa)/3m. The total number of 
atoms N = J|n(r)| 2 dr and the total magnetization 
M = J [|n + (r)| 2 — |n_(r)| 2 ] dr are conserved quantities. 
The Zeeman-energy shifts for each component can be cal- 
culated using the Breit-Rabi formula [l5| 
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where -Ehfs is the hyperfine energy splitting at zero mag- 
netic field, a — (gj + gj)/j,BB/ -Ehfs, where p,s is the 
Bohr magneton, gj and gj are the gyromagnetic ratios 
of electron and nucleus. The linear part of the Zccman 
effect does not affect the condensate evolution, except 
for a change in the relative phases [l6| and so we re- 
move it with the transformation \P ± — > 4'±cxp(— iE±t), 
*o -* * exp[-i(E + + E_)t/2]. We thus consider only 
the effects of the quadratic Zeeman shift, SE = (E + + 
£L — 2E )/2 fa a 2 Enps/16, which is always positive. 



III. MODULATIONAL INSTABILITY OF 
HOMOGENOUS STATIONARY STATES 

First, we are interested in the stability analysis of the 
homogeneous condensate and consider the case of van- 
ishing potential, V(r) = 0. We look for the homogeneous 
solutions in the form ipj = -/nje l/l3 ' i+ * . The "phase 
matching condition" for Eqs. ([T]) gives /i+ + \i- = 2/xq- 
We find that both in the case of B = and in the case 
of M = 0, the steady state fulfills the stronger condition 
//+ = H- = Ho- However, if both magnetic field and mag- 
netization are nonzero, which is the case in real experi- 
ments, the chemical potentials will be different, satisfying 
the less stringent phase matching condition. 

We define the density fraction in each component as 
Pj = rij/n. If we assume that all three spin components 
Pj are nonvanishing, the relative phase between them, 
6 = 26o — + — 6- , can take the value or ir only. We will 
describe the corresponding stationary states as phase- 
matched (9 — 0) and anti-phase-matched {6 = tt). Note 
that both types of states can in general exist in both fer- 
romagnetic and polar condensates [l?], EH • However, 
phase-matched states are energetically favorable in fer- 
romagnetic condensates, and anti-phase-matched states 
in polar condensates [3, [IB]. For that reason, they 




FIG. 1: (Color online) Diagram of existence of station- 
ary states in spin-1 condensates. In addition to the three- 
component phase-matched and anti-phase-matched solutions 
shown on the diagram, two-component solutions with po = 
exist for arbitrary M, and one-component solutions pj = 1 
exist with j = — ,0, +. Dark-shading (top) - anti-phase- 
matched state, light-shading (bottom) - phase-matched state, 
cross-hatching - both states. The dotted and dashed lines at 
M — indicate the absence of a phase-matched or anti-phase- 
matched state respectively. 



were named ferromagnetic and polar states respectively 
in Ref. [l6| . In Fig. []] we present the existence diagram 
for three-component homogeneous stationary states. For 
generality we also include the results for ferromagnetic 
condensates. The ferromagnetic condensates, such as 
87 Rb, occur in the lower half (where c 2 is negative), while 
polar condensates, such as 23 Na, occur in the upper half 
(c 2 is positive). There is clearly a region of coexistence of 
anti-phase-matched and phase-matched states for a po- 
lar condensate in nonzero magnetic field. In addition, a 
two-component solution with p = 0, and one-component 
solutions with pj = 1 exist. Our results are in agreement 
with the previous analysis of homogeneous ground states 

B 

The energy density is related to the Hamiltonian of the 
system, from which Eqs. (p} are derived, by H = J Edr. 
In addition to the anti-phase-matched ground state [3| , 
the polar condensate in the coexistence region of Fig. [T] 
has an excited phase-matched state corresponding to the 
energy maximum at 9 = 0. This state is stable with re- 
spect to spatially homogeneous spin mixing, because the 
possible dynamical trajectories in the (po,0) plane cor- 
respond to a constant energy value, hence both minima 
and maxima are stable. 

The stability properties of these states change when 
we consider the possibility of a spatial, or modulational, 
instability (MI) . We calculate the growth rate of the Bo- 
goliubov modes jl9j . 
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After substituting the above to Eq. ([T]) we obtain a set 
of equations for the vector z = uq, u_, v*L, Vq, v*_) 
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TABLE I: Stability of spin-1 condensate states in absence and 
presence of the magnetic field: PM - phase-matched, APM - 
anti-phase-matched, X - state does not exist. 1 A family of 
stationary states. 2 Neutral stability with respect to spatially 
homogeneous spin mixing. 



in the form dz/di = iAz, where A is a 6th-rank matrix 
(20| . For an equilibrium state, it is possible to eliminate 
pij and SE from A, expressing it in terms of rij, cq, C2 
and k only for specific 9. The Bogoliubov modes are the 
solutions of the characteristic equation det(A — Tilo) = 0, 
with ijj being the eigenfrequency of the excitation. The 
form of this equation is too cumbersome to present here, 
in contrast to the B — case considered before [ll|. In 
general, one has to use numerical methods to obtain a 
set of solutions. 

The general numerical results are presented in a sys- 
tematic way in Table fl] These results have been calcu- 
lated for the specific cases of the 87 Rb condensate (upper 
half) and 23 Na condensate (lower half) with scattering 
lengths given in [l6| and hyperfine energy splitting given 
in [2l| . As we are considering the homogeneous case the 
result is applicable to one-, two- and three-dimensional 
condensates. We do not consider here the trivial case of 
p± = 1 states, which have extremal value of magnetiza- 
tion M = ±N. All the cases presented in this table were 
double checked by direct numerical integration of Eq ^ 
(see Sec. II VI) . The results cover the area of Fig. [1] and do 
not apply to the case of very high magnetic field, when 
the quadratic Zeeman energy dominates. 

We find that our results are in agreement with the 
existing data for the vanishing magnetic field case. In 
Ref. the authors found that all the condensate equi- 
librium states are stable, with the exception of one state 
in the ferromagnetic condensate, which corresponds to 
"ferro- APM", M = state in Table H In Ref. 0, it 
was also found that ground states of both ferromagnetic 
and polar condensates ("ferro-PM" and "polar- APM'' ) 
are modulationally stable. The authors of Ref. [ll[ were 
considering mainly non-equilibrium (spin- mixing) states, 
but their general conclusion was that all polar condensate 
states are dynamically stable, and ferromagnetic are not. 
Here we report that in fact specific ferromagnetic con- 
densate states are stable, while polar condensates may 
become unstable. Furthermore, as discussed below, we 
see that dynamic domain formation may occur for polar 
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FIG. 2: (Color online) Modulational instability growth rate 
in a phase-matched steady state of a sodium condensate (cir- 
cles), versus the quadratic Zeeman shift in the case M = 0. 
The continuous line is a square fit to the numerical data. The 
upper shaded area shows the range of k vectors corresponding 
to the unstable (imaginary) frequencies. 



condensates, and lead to convergence to new stationary 
states. 

As one can see in Table HI the magnetic field af- 
fects stability of polar condensates. We investigate this 
phenomenon in detail by calculating the growth rate 
K = Im(w) corresponding to unstable Bogoliubov modes. 
In Fig. [2l we present results for a "polar-PM" state of 
a sodium condensate as a function of the magnetic field 
strength. The growth rate is proportional to the square 
of the quadratic Zeeman shift, which is in turn propor- 
tional to the square of the magnetic field strength. Hence, 
one has to apply a relatively strong magnetic field to ob- 
serve MI on a reasonable time scale. Another interesting 
feature is the range of wavevectors k corresponding to 
unstable modes. In contrast to the typical case where 
this range starts from k = 0, here the unstable region 
begins at a nonzero minimum value. This type of "op- 
tical mode" branch has been reported before in the case 
of parametric optical solitons [2Q| • 



IV. DYNAMICS OF THE CONDENSATE IN A 
CIGAR-SHAPED TRAP 

We now consider the implications of these results 
in the experimentally relevant case of a 23 Na conden- 
sate localized in a cigar-shaped harmonic trap, V(r) — 
^mu}^_(y 2 + z 2 ) + jmuj 2 x 2 . Specifically we consider the 
case luj_ » uu in which the Fermi radius of the trans- 
verse trapping potential is smaller than the spin healing 
length, and nonlinear energy scale is much smaller than 
the transverse trap energy scale, which allows us to re- 
duce the problem to one spatial dimension fl6l . |22| . Fol- 
lowing standard dimensionality reduction procedure we 
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obtain the one-dimensional model 
dip± 



(a) PM 



(b) APM 



ih 



dt 
dt 



C + c 2 {n± + n Q - n T ) 
£ - SE + c 2 (n+ + n,-) 



i/j + 2c 2 ip + ip-ipQ 



(4) 



where C = —(ft, /2m)d /dx +cqu+ ^muj^x and the in- 
teraction coefficients have been rescaled and now include 
the transverse trap frequency, Co = / ihujj_(2a2+ao)/3 and 
c 2 = AHuj±(a2 - Oq)/3. 

The experimental scenario we consider here consists of 
several phases. Initially, the condensate is prepared in 
the m = — 1 ground state [12]. Next, short microwave 
field pulses are applied to transfer the atomic population 
to the desired state [23| . We consider two phase- 
matched state (6 = 0) with p+,a- = 0.351,0.3,0.349, 
and a anti-phase-matched state (9 = n) with p+,o,- = 
0.4,0.01,0.59. Simultaneously, the magnetic field is set 
to the value of 175 mG [12]. The results of the corre- 
sponding numerical simulations of Eqs. ^ are presented 
in Fig. [3] The MI develops after tens of milliseconds, 
and leads rapidly to spin-domain formation. We see what 
appears to be initial oscillations, followed by instability 
dynamics leading to an apparent oscillating state. In 
the case of the phase-matched initial condition and al- 
most zero magnetization we ultimately see conversion to 
the "00 component, as observed in experiment fl3j |. In 
the anti-phase-matched initial condition however, with 
nonzero magnetization, the spin domains become well- 
defined and persist in the dynamics. We have verified 
that the nonlinear energy at peak density Con max /2 is 
much smaller than the transverse energy separation ftuj±, 
which justifies the use of reduced ID GPE (22j. 
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FIG. 3: (Color online) Spin-domain formation in a 23 Na con- 
densate confined in an optical harmonic trap. Top panels: 
Evolution of no (darker shading higher density) . Bottom pan- 
els: densities at given times. (a,b,c) phase-matched initial 
state, p+,o,- = 0.351,0.3,0.349; (b) anti-phase-matched ini- 
tial state, p+,o,- = 0.4, 0.01, 0.59. The thin dotted line corre- 
sponds to the total condensate density n, and the dashed, 
solid and dotted lines correspond to n+, no and n_ re- 
spectively. Parameters are N = 3.7 x 10 3 , B = 175 mG, 
lo± = 2tt x 10 3 Hz, wii = 2tt x 32 Hz. 



V. STATIONARY SPINOR STATES 

The existence of MI suggests that domain-type station- 
ary states may be expected in the trap, just as plane wave 
instability and solitons are typically found together in op- 
tics. Indeed, as can be seen in Fig. [4] we find that in the 
stationary picture the profiles always break the single- 
mode approximation for both the phase-matched and 
anti-phase-matched states (as found in Ref. [l4| for the 
case of a anti-phase-matched state in a polar condensate). 
As predicted by the MI analysis, an initially smooth pro- 
file will therefore become modulated with the ensuing in- 
stability dynamics reflecting the nearby stationary state 
profiles. For instance comparing Fig. [3^b) and Fig. [4ja) 
we see that profiles similar to the stationary states ap- 
pear in the evolution. Significantly, but unsurprisingly in 
a polar condensate, we find that the phase-matched state 
is generally unstable, while the anti-phase-matched state 
is stable. The stability of these states appears to reflect 
the ultimate dynamics of the condensate with the phase- 
matched state breaking up while the anti-phase-matched 
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FIG. 4: (Color online) Examples of (a) phase-matched and 
(b) anti-phase- matched stationary states in a 23 Na conden- 
sate. Shown are the total density n (thin dotted line) and 
component densities n+.o,- (dashed, solid and dotted lines 
respectively) . 



state has highly persistent spin domains. A complete 
analysis of the stability of the stationary states will be 
presented elsewhere. 
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VI. CONCLUSIONS 

We have demonstrated that an antiferromagnetic spin- 
1 condensate can undergo a novel type of spatial modula- 
tional instability followed by subsequent spin-domain for- 
mation in the presence of a homogeneous magnetic field. 
We have employed realistic conditions to demonstrate, 
with the help of numerical simulations, that this novel 
modulational instability can be observed in a sodium con- 
densate confined in an optical trap potential and that the 



ensuing instability dynamics connect with the stationary 
states in the trap. 
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